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Цель данного исследования — определение условий 
повышения эффективности функционирования балан- 
сировочного оборудования. В качестве задачи рассмат- 
ривается повышение точности и снижение трудоемко- 
сти балансировки ротора. В основе исследования — 
модель ротора на балансировочном станке, рассмотрен- 
ная с точки зрения теории колебаний как динамическая 
система с двумя степенями свободы. В результате про- 
веденных изысканий, во-первых, обоснована необходи- 
мость использования в модели ротора векторных значе- 
ний коэффициентов влияния дисбаланса в каждой из 
плоскостей коррекции на колебания опор ротора на ба- 
лансировочном станке. Во-вторых, представлена мето- 
дика их экспериментального определения. Полученные 
результаты могут быть использованы при создании ба- 
лансировочного оборудования (балансировочные стан- 
ки и балансировочные комплекты) нового поколения. 
Проведенное исследование представляет собой разви- 
тие современной теории и практики балансировки, 
направленное на повышение адекватности динамиче- 
ской модели ротора. 
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адекватность, векторные коэффициенты влияния. 


ТБе \огк оБ]есйуе 1$ ю д&егтше Фе сопа1юоп$ Гог 1т- 
ргоуш> фе орегайп> ейслепсу оГ Фе Баапсте еда1р- 
тепе. ТБе ргоет 15$ ю шсгеазе ассигасу ап ю гедисе 
]1абог пиепзКу ог Ше гоог Ба!апсшо. А гоог тоде| оп Ше 
Ба]апсте тасфше сопз1Аеге4 ассогте {10 Фе угаНоп 
Феогу аз а Чупаптис зузет \/ИЫ (мо Дестеез ог Неедот 
Гогил$ Фе Баз1$ оЁ Ше гезеагсВ. Аз а гези о? Фе шуезИза- 
(оп$, Нгзё, фе песеззИу оЁ аррушх уесюг шЙчепсе сое{- 
ЙслепЕ уашез ог Фе гоюг то4е| ипра]апсе ш еасВ оР Фе 
соггесйоп р]апез оп Фе угайоп$ оЁ Фе гоюг ро[ез$ оп Фе 
Ба!апсег 1$ ргоуе4. Зесоп у, Фе (есбтауие оЁ Фе ехрег- 
ппеша| деегиитайоп 1$ ргезещед. ТБе геза $ оМашеа 
сап Бе изе4 {о сгеме Байапстпе едитртепе (Баапсте та- 
сыпез апа Баапсте $е($) оЁ пе\ сепегайоп. ТБе зшау 15$ 
Фе деуеортепЕ о? Фе то4егп Феогу апа ргасисе Баапс- 
12 аппед а паргоуше Фе адедаасу оЁ Фе дупапис гоог 


подЕ]. 


Кеумога$: Ба!апсше едилтртепь еЁ@слепсу, Дупаписз, 
Гопсйопа| п1о4е], адеаиасу, уесюг шЙчепсе сое слеп. 


Введение. В современной производственной практике приходится решать такие проблемы, как снижение 
вибрации и шума машин, повышение их надежности и качества функционирования, улучшение условий труда персо- 
нала, управляющего работой машины. Эффективным средством решения подобных задач является балансировка ро- 
торов — в том числе и высокооборотных карданных передач автомобилей. Этим объясняются высокие требования, 
предъявляемые к балансировке роторов. В первую очередь речь идет о качестве ее исполнения (на эксплуатационной 
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частоте вращения остаточные дисбалансы ротора в плоскостях его коррекции не должны превосходить допустимые 
значения, регламентируемые технической документацией). Кроме того, следует учитывать трудозатраты балансиров- 
ки в производстве машин и при их ремонте (минимизация числа балансировочных циклов как основа снижения трудо- 
затрат балансировки). 

Реализация этих требований напрямую зависит от адекватности механико-математической модели балансиру- 
емого изделия. Использование такой модели на балансировочном станке позволяет установить связь между характе- 
ристиками колебаний опор балансируемого изделия и дисбалансами в плоскостях его коррекции с учетом их взаимно- 
го влияния на колебания опор. 

В первых станках для динамической балансировки роторов (например, станки рамного типа конструкции 
Б. В. Шитикова [1-3]) задача решалась механическим исключением взаимного влияния плоскостей коррекции. С этой 
целью последовательно включалась и отключалась податливость шарнирных опор рамы, совмещенных с плоскостями 
коррекции балансируемого изделия. При этом использовалась простейшая механико-математическая модель баланси- 
руемого изделия на станке 


А; =К;хО;; 1=12, (1) 
где А; = А; Хе"?! — амплитуда А; и фаза ф; колебания податливой 1-й опоры рамы станка при абсолютно жесткой 
другой ее опоре; К; — скалярный коэффициент влияния, определяемый экспериментально; 0; — оценка дисбаланса в 


1-Й плоскости коррекции, сонаправленного с А; (р; = р; хе, где О; — значение, а ф; — угол дисбаланса р;). 

Такая модель имеет низкую адекватность и требует большого количества балансировочных циклов, что объ- 
ясняется необходимостью последовательных замеров и корректировок дисбалансов отдельно в каждой из плоскостей 
коррекции. 

В функционировании появившихся позже станков маятникового типа использовалась модель [4]: 


А} =[ К, хр; 1=Ь2;=12, (2) 


где [АП —щ матрица-столбец амплитуд и фаз синхронных колебаний опор 1 = 1, 2 ротора на станке; {р,} — матрица- 
столбец искомых дисбалансов в плоскостях коррекции ] = 1, 2 отстоящих от левой и правой опор ротора на расстоя- 


ния 1 и С соответственно; [К; й — матрица скалярных коэффициентов влияния: 


[-А А 


С 1-С7 ’ 


К, |=(1-А-С)х 


где [ — расстояние между опорами ротора на балансировочном станке. 
Заметим, что здесь и ниже обозначения матриц столбцов раскрывают и по [5] как {Х;} = {Х, Х.....Х,....Хи} = И Х.... 
Х....Х,) при [= 1...п. 

В отличие от (1), модель (2) позволяет получить оценки дисбалансов {р,} сразу в двух плоскостях коррекции 
за один балансировочный цикл. Однако эти оценки имеют низкую достоверность, т. к. матрица (3) предусматривает 
статическое распределение дисбалансов между опорами, не учитывая динамические свойства колебательной системы 
ротора на балансировочном станке. Поэтому, сокращая число балансировочных циклов определения и корректировки 
{2;}, модель (2), как и модель (1), не имеет достаточно высокой адекватности. 

Повышение эффективности функционирования станков маятникового типа привело к необходимости исполь- 
зования механико-математической модели (2) с матрицей скалярных коэффициентов влияния [4, 6-8]: 


Че Сы 
К око! (4) 
21 р, 
Элементы которой определяются тарировкой измерительной системы станка на частоте вращения, принятой для ба- 
лансировки изделия. 


Если речь идет об универсальных балансировочных станках, такая матрица строится для заданных совместно 
диапазонов масс балансировочных изделий и частот их вращения при балансировке. 
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Модель (2) с матрицей (4) коэффициентов влияния лишь количественно учитывает динамические свойства 
колебательной системы балансируемого изделия на балансировочном станке. Поэтому дальнейшее повышение адек- 
ватности механико-математической модели (2) целесообразно искать с использованием матрицы векторных коэффи- 
циентов влияния 

к, | = Ки Кр | (5) 
Ку Кр 
определяемых экспериментально для принятой частоты вращения ротора при его балансировке. Это предложено в 
работе [9], цитируемой по [10]. Однако на практике данный подход не получил развития из-за отсутствия теоретиче- 
ских обоснований возможности и необходимости его использования, а также из-за сложности его реализации с помо- 
щью предложенного в указанных работах графоаналитического метода. 


Ниже представлены необходимые теоретические обоснования и практически реализуемая методика опреде- 


ления векторных коэффициентов влияния Ге. [=1,2; /=1,2, обеспечивающих наибольшую адекватность ис- 


пользования механико-математической модели динамики неуравновешенного двухопорного ротора на балансировоч- 
ном станке: 


А} = К, жи =12; 1 =1,2. (6) 
Теоретический анализ коэффициентов влияния. На рис. | представлена динамическая модель идеального 
однородного жесткого двухопорного ротора с податливыми и демпфированными опорами на балансировочном станке. 


Ах 





Рис. 1. Модель для анализа колебаний опор ротора на балансировочном станке 


от проявления дисбаланса Д в левой плоскости коррекции (] = 1) 


Модель имеет две степени свободы, описывающие вынужденные колебания вращающегося со скоростью © 


ротора в горизонтальной плоскости хс. Эти колебания генерируются усилием Вх” от дисбаланса р=Шхе*%? со 
значением ШО. Здесь фр — угол идеально уравновешенного ротора, установленного в плоскости коррекции, отстоящей 
от левой опоры на расстояние А. На рис. | показаны положения колебательной системы: статически равновесное 
(сплошная линия) и текущее (пунктирная линия). Здесь также обозначены: 
— [ — как и ранее, расстояние между опорами ротора; 
— 5 — положение центра масс ротора; 
— К.,К> — жесткости 1-й и 2-Й опор ротора на балансировочном станке; 
— с1, С2 — коэффициенты диссипативности опор; 
— х., Х2 — перемещения опор, которые принимаем за обобщенные координаты. 

Для составления уравнений движения в обобщенных координатах используем уравнения Лагранжа второго 
рода с учетом диссипативных сил [11]: 


а[от | ое о(Т-И) 


ВВ ОО ЕЕ + = 


в которых Т, И — выражения для расчета кинетической и потенциальной энергии; ЁР — диссипативная функция; 
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0; — обобщенные силы. 


Выражая через обобщенные координаты х,,х› перемещение х центра 5 масс ротора и поворот ф его оси как 


ИЕ 
ОЕ, 
2 | 
раскрываем входящие в (7) выражения: 


Т-тх +51 ф тата х(и-№): 


1 2 1 2 1 .2 1 в. 2. 
= К хх +- Кхм: Е =-ахх +-с хм»: 
мм То ^ 72 ри о 62 >72 





О =Рхо”х хс0$ (@0#+Фр }; 


(9) = хо? к* хо + Фр), (8) 


в которых ит — масса ротора; / — момент инерции ротора относительно оси у. 
Выполнив дифференцирование выражений (8) по форме уравнений Лагранжа (7) и используя для решения 
полученных уравнений метод комплексных амплитуд [11], представляем эти уравнения в матричном виде: 


Саре ера = (Ве, Й 


1 Г 1 1 
м т 
где [А] =: Е = — матрица инерционных коэффициентов; [В] = фаз [с со | — диагональная матрица 


диссипативных коэффициентов; |С | = Фа [© К | —щ диагональная матрица квазиупругих коэффициентов; 
{4} = {. Х.}; {9} = {хх}; {4} = (х, х2} — матрицы-столбцы обобщенных ускорений, скоростей и перемещений 


2х ©” 
соответственно; Е} = 





х }(1 — А) А] — матрица-столбец обобщенных сил. 


При решении уравнения (9) в виде {4} ={}}х е'®' подстановкой входящих в него величин получили выраже- 
ние для расчета матрицы комплексных амплитуд в виде 
-1 
2 . 
{7} = ([с]-® х[А]+5х[ В] Е, (10) 
где степень —1 указывает на обращение матрицы, заключенной в круглые скобки. 


Раскрывая комплексные амплитуды через их действительные (а;) и мнимые (Ъ;) слагаемые и используя фор- 
мулу Эйлера, можно записать: 


{1} = {а +в а, +, } = {А же А; же“ |, (11) 
где А; = \ а? + ; а; = 81510 (Ь, ) агссо$ т амплитуда и фаза колебаний 1-й опоры ротора на балансировочном стан- 
ке (1 = 1,2). 

Элементы матрицы {}} по (11), определенные по (10), — комплексные величины, т. к. элементы квадратной 
матрицы 
ы —1 
[Е] = ([<]-в? х[А]+х[ В], (12) 


описывающей собственные свойства колебательной системы ротора на балансировочном станке, являются комплекс- 
ными числами из-за учета объективно существующего демпфирования в опорных элементах ротора на станке. 
Раскрытием входящих в (12) матриц получаем 
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—1 
К -®х ее + ХС! —0°х и 
т 4 р 4 [2 
[Ш] = 
2 (1 1 2 [1 Г). 
—0х| —т-— К -®х| —т+— |Н1юхс, 
4 р 4 [2 
что после обращения матрицы позволяет найти: 
ый [и ое, 0х а 
1 4 РЁ 4 Г 
[2] = | (13) 
ее | Г. > | > 2 | | 
Фо х| —т-—— К -®х| —-т+— | +юхс, 
4 Г Г 
где 
2 
Че. = К -ю®’х ы юж 1х4 К — 6х о юж ®'х ЕЯ (14) 
1 4 12 1 2 4 12 я, 4 12 


— определитель матрицы [[]. 
Если в рассматриваемой на рис. 1 плоскости коррекции (ее обозначили как } = 1), расположенной на расстоя- 


нии А от левой опоры ротора, установлен единичный дисбаланс О, имеющий значение Д = 1 и угол фр = 0, то опреде- 
ленные по (10) комплексные амплитуды будут представлять собой искомые векторные коэффициенты влияния дисба- 
ланса в ] = 1-й плоскости коррекции на колебания опор 1 = 1,2 ротора: 


й. 
(Ки кы| (ки хе К. ке} = "Е х{(1-А)А}. (15) 


На рис. 2 представлена динамическая модель того же, что и на рис. 1, ротора, имеющего те же собственные 


и ь и В —1 
свойства колебательной системы на балансировочном станке (они описываются той же матрицей [2] по (12)). 


Ах 





Рис. 2. Модель для анализа колебаний опор ротора на балансировочном станке 


от проявления дисбаланса О в правой плоскости коррекции (] = 2) 


Колебания опор этого идеально уравновешенного ротора при его вращении на балансировочном станке воз- 


буждаются лишь дисбалансом О, установленным в плоскости коррекции /=2, расположенной на расстоянии С от 
правой опоры 1 = 2 ротора. 

Колебательный процесс модели на рис. 2 будет описываться тем же уравнением (9), в котором матрица — 
столбец обобщенных сил должна представляться как 


{в} = 


Принимая О =1, Фр=0 и подставляя (16) в (10), получаем уравнения для расчета коэффициентов влияния 


) 2 
т ео] (16) 


дисбаланса в плоскости коррекции ] = 2 на колебания опор 1 = 1,2 ротора на балансировочном станке: 
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2 
Е — (к хе®> К, хе®» и -х[Е х1С(1-С), (17) 


в котором [1] 'та же, что ив (15). 

Анализ теоретических зависимостей (15), (17) и входящих в них соотношений (13), (14) позволяет утверждать 
следующее: 
— учитывая наличие демпфирования колебаний опор ротора на балансировочном станке, коэффициенты влияния 
дисбалансов в плоскостях коррекции /] = 1,2 на колебания опор 1 = 1,2 ротора необходимо рассматривать как вектор- 


ные величины (К; =К, же 1 =1,2 ‚1=1,2 ); 


— значения К, и углы В; коэффициентов влияния К‚„ существенно зависят от скорости ® балансировки ротора 


на станке. От в зависят также характеристики [А колебаний опор неуравновешенного ротора на балансировочном 


станке. Поэтому для использования модели (6) при балансировке ротора элементы матриц [А И К. Должны 


быть определены на одной и той же скорости @® вращения ротора. При этом значение в балансировки жесткого ротора 


не влияет на определяемые посредством (6) характеристики {р _ его неуравновешенности; 


— векторы К; практически не совпадают как по значениям, так и по углам для любых сочетаний [и /. Лишь теорети- 


чески не исключен случай, когда К =К,, с! =с> , А = С ‚, при котором имеет место а = К Ко — К» : 
— значения и углы векторных коэффициентов влияния для принятой скорости © балансировки заданной конструкции 
ротора на определенном балансировочном станке зависят только от координаты, определяющей положение плоскости 
коррекции относительно одной из опор ротора. 

Поясняя последнее, отметим, что при подстановке в (15) задаваемого значения 4, определяющего положение 





плоскости коррекции / = | относительно опоры 1 = 1, рассчитывается матрица-столбец га К. коэффициентов вли- 


яния. Если в том же уравнении (15) принять А = [ — С — координату плоскости коррекции / = 2 относительно 1 = 1-й 


опоры, то результатом расчета будет матрица-столбец {К К». коэффициентов влияния. Те же результаты получа- 


ются и с использованием уравнения (17). В этом случае в расчет вводятся значения С — расстояния от плоскости кор- 


рекции / = 2 до опоры 1 = 2 ротора на станке, и таким образом находится матрица-столбец {К К, коэффициентов 


влияния. Принимая в том же уравнении С =[-А — расстояние от опоры 1 = 2 до плоскости коррекции ] = 1, нахо- 
дим матрицу-столбец га Ко коэффициентов влияния. 


Таким образом, по каждому из соотношений (15), (17) могут быть найдены коэффициенты влияния дисбалан- 
са в любой плоскости коррекции (не только в ] = 1, 2) на колебания опор {= 1, 2 ротора на балансировочном станке. 
Для этого задается либо расстояние А (рис. 1 с использованием уравнения (15)), либо расстояние С (рис. 2 с использо- 
ванием уравнения (17)). 

Методика определения коэффициентов влияния. При балансировке конкретного изделия находить вектор- 
ные коэффициенты влияния аналитическим методом нерационально из-за трудоемкости и неточности определения 
величин, входящих в (15) или (17), и, как следствие, значительной погрешности результатов расчета. Аналогичным 
образом объясняется отказ от использования векторных коэффициентов влияния с графическим методом их опреде- 
ления [9, 10], о чем упоминалось выше. 

Излагаемая ниже методика натурного определения векторных коэффициентов влияния для жесткого дву- 
хопорного ротора строится на ряде допущений. 


— Амплитуды А;, фазы @; колебаний А; =А,; хе” опор 1=1,2 ротора и частота его вращения при балансировке 


определяются только характеристиками неуравновешенности ротора. 
— Неуравновешенность ротора полностью определяется искомыми и корректируемыми при балансировке дисбалан- 


Ф. 


—— 7 . 
сами О = 72) ‚хе в двух заданных плоскостях его коррекции (] = 1,2). 


— При анализе малых колебаний опор ротора от различных источников используется принцип суперпозиции. 
Возможность принятия этих допущений при моделировании неуравновешенности и колебаний жестких роторов на 
балансировочном станке доказана существующей практикой балансировки. 
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Рассмотрим ротор, имеющий в заданных плоскостях коррекции ]/=1,2 искомые дисбалансы р, =Бх е'®' 
р, =, хе“. С приведением этого ротора во вращение на балансировочном станке колебания А, его опор будут 
описываться выражением (6), раскрываемом как 
А=КихВ+Крхр» ; 
А =КухЬ +К»хр, . (18) 
Установим в плоскости коррекции / = 1 пробный дисбаланс р =Рих е®" с заданными значением Ри: и уг- 
лом фп/. После приведения во вращение колебания опор этого ротора будут описываться: 
АА =Ких(В Ри )+К хр» ; 
А, =Кых(Р+Р}+К» хр. (19) 
Вычитая из (19) выражения (18), находим 
А! -А = Ау хе — Ал хе“ =К. хр же“; 
Аз! — Аз = А хе" — А, же =К хр, хе. 
Отсюда 
А! хе А, хе 


О ” 
Ри хе й 


Г. 1 — 
Ки =Кухе"" = 


1421 __ 102 
А хе А хе 


К =К.. хейл = 
21 р х Фит 
п] е 


(20) 


Удалим из плоскости коррекции ] =1 дисбаланс )‚, и установим в плоскости коррекции ] =2 пробный дисбаланс 


Р.› =О,› хе? с заданными значением О, и углом Ффп2. Опишем колебания опор этого ротора на станке: 


А =КихР +К»›х(Р, +, } (21) 


Вычитая из (21) уравнения (18), окончательно получаем: 
12 10 
А} хе Ау хе '. 


к — Вр? — 
Ко =Крхе ® = = | 
О, хе 


1952 _ 10,5 
А5> хе А хе 


к _— В>2 _ 
К =Кохе® = т. 
О, хе 


(22) 


Выражения (20), (22) определяют векторные коэффициенты влияния через заданные характеристики пробных 


дисбалансов О, О, и определенные замерами на балансировочном станке характеристики (амплитуды и фазы) ко- 


лебаний опор 1 = 1, 2 ротора с пробными дисбалансами и без них. 


Формируя по результатам расчетов матрицу К. по (5), находим дисбалансы р, р, балансируемого ротора в при- 


нятых плоскостях его коррекции: 


РА=[Е, | {А}. (23) 


При необходимости балансировки партии изделий одного типоразмера (как это требуется, например, при ба- 
лансировке карданных передач в производстве и ремонте автомобилей одной модели) матрица [К В построенная по 
изложенной методике для одного изделия, «обучает» измерительную систему станка для использования ее при балан- 
сировке всех изделий партии. 

В лаборатории «Динамика приводов, вибрации и диагностика качества машин» Донского государственного 


технического университета был поставлен эксперимент. В ходе работы была реализована балансировка двухшарнир- 
ной карданной передачи 2121-2201012-04 автомобиля «Нива» по существующей технологии с использованием матри- 
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цы скалярных коэффициентов влияния (4) и балансировка той же передачи по новой технологии с использованием 
матрицы векторных коэффициентов влияния (5). С учетом контрольного измерения первая технология была реализо- 
вана за 5 балансировочных циклов, вторая — за два цикла при одних и тех же значениях допустимых дисбалансов. 

Заключение. Для разработки балансировочного оборудования нового поколения и прогрессивных техноло- 
гий балансировки следует использовать механико-математические модели. В данном случае речь идет о моделях ди- 
намики ротора с векторными коэффициентами влияния дисбалансов в плоскостях коррекции на характеристики коле- 
баний опор. Целесообразность такого подхода теоретически обоснована. Предложенная методика определения век- 
торных коэффициентов влияния обеспечивает получение достоверных результатов и легко реализуется на компьюте- 
ре. Обеспечивая более высокую адекватность, модель с векторными коэффициентами влияния позволяет повысить 
точность определения дисбалансов балансируемого изделия в плоскостях коррекции, уменьшая тем самым необходи- 
мое число балансировочных циклов и трудозатраты балансировки. 
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